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Ih this paper the influence of the eddy currents on a telephone relay of the 
electromagnetic type on closing is investigated. The same relay as in Parts 1, 2, 
and 3 of the study is considered, i.e., an ’*L-armature” relay with a round core 
and a rectangular yoke, where the magnetic leakage is taken care of by the 
simplified field picture introduced in Part 1. The relay is assumed to have eddy 
currents in core and yoke. The armature is kept fixed and the magnetic reluctances 
are supposed to be constant. 

The concept of an equivalent winding, representing the eddy current paths in a 
magnetic circuit, is subjected to a general theoretical treatment. It is shown that 
an equivalent winding should give useful results on closing but cannot be employed 
on breaking. 

For the relay under consideration the equivalent winding parameters on closing 
are deduced. Numerical results are reported. The importance of the magnetic 


skin effect is pointed out. 
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CHAPTER 1 


Introduction and summary 


In the present Part 4 of ““A Study of Telephone Relays” the same type of relay is con- 
sidered as in the earlier Parts 1, 2, and 32-4, i.e., a relay of the L-armature type with a round 
core and a rectangular yoke. As before, the armature is kept fixed, the magnetic reluctances 
are assumed to be constant, and the magnetic leakage is represented by a leakage flux, 
running perpendicularly between core and yoke. 

The present paper and Part 3 are closely connected inasmuch as both papers discuss the 
same aspect of relay theory, viz., the influence of the eddy currents arising in the metallic 
parts of the magnetic circuit. As in Part 3 we assume eddy currents in both core and yoke. For 
some general points of view reference is made to Part 3, Chapter 1. 

Employing the concept of operational reluctances, we could in Part 3 derive without 
difficulty general operational expressions for the magnetic core flux on closing or breaking, 
and for the magnetizing primary current on closing. On breaking—the problem studied in 
detail in Part 3—it was then possible to deduce from the operator of the magnetic flux (only 
the armature end flux was: considered) two different expressions for the corresponding time 
function: Heaviside’s expansion theorem gave us an exponential function series, suitable for 
large values of the time /;’a series expansion of the operator in powers of p-1/4 gave us a 
series in powers of 1/4, suitable for small rt. These two solutions were used for studying 
quantitatively an actual relay. The results yielded by the study were in excellent agreement 
with actual oscillograph records of the flux decay. 


In this paper we turn to a detailed study of the problem of closing, trying to obtain suitable 
expressions for the primary current and the magnetic flux. Unfortunately, in this case the 
method of Part 3 seems to be excluded on account of numerical complications (as to the 
expansion theorem solution see the remark following eq. (119) in Part 3; as to the power 
series solution, given by. eq. (124), this solution can only be used for extremely small 
values of t, and further terms of the series are difficult to obtain, as they require a better 
approximation of the quantity 72than the one given by (72)). 

A new approach is necessary. In the following it will be shown that on closing sufficiently 
accurate solutions can be obtained by the method of replacing the eddy current paths by an 
equivalent winding. 


The idea of an equivalent winding was mentioned already in Part 3, Chapter 1. It was 
pointed out that some earlier writers’:5»® have made use of this idea, interpreting more or less 
explicitly the equivalent winding, as a physical wire winding, wound around the relay core 
(cf. also a paper by H. Hess and E. WeserR’®). 
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The results of Part 3 strongly supported the assertion that the breaking problem cannot be 
studied by the aid of a single equivalent winding. This procedure leads to an exponential 
decrease of the flux-time curve, an outcome which is very far from the actual course of events. 
Other and more suitable methods could, however, be found. 

For the closing problem, on the other hand, where—as explained above—these methods 
break down, the concept of an equivalent winding is useful. One of the aims of the present 
paper is to develop a satisfactory general theory of the equivalent winding. 


Chapter 2 deals with the basic principles of the equivalent winding, and gives general 
expressions for its parameter values, as depending upon the eddy current distribution. In 
Ckapter 3 these expressions are used for an approximate calculation of the equivalent winding 
parameters of a relay on closing. 

In Chapter 4 the results are tested on a case which can also be exactly solved, viz., a relay 
without magnetic leakage and with eddy currents in the core only. The comparison with the 
exact solution turns out satisfactorily. In Chapter 5 the equivalent winding is then used for a 
quantitative study of the closing of a relay with magnetic leakage and with eddy currents in 
core and yoke. In order to facilitate the computations, a list is finally given of the Principal 


symbols and formule. 


CHAPTER 2 


General theory of the eddy current equivalent winding 


This chapter gives a general theory of the equivalent winding of the eddy currents in a 
magnetic circuit. The theory is presented as an application to our telephone relay. It is quite 
general, however, and can in principle be applied to any magnetic circuit. 

We consider a relay according to Fig. J] on the next page (cf. Part 1, Fig. 4; Part 2, Fig. 9; 
Part 3, Fig. 2) with a cylindrical core and a rectangular yoke. The core and the yoke are 
assumed to be the only electrically conducting parts of the magnetic circuit. A primary 
magnetizing wire winding with N, turns is evenly distributed over the whole length of the core. 
The relay has a magnetic leakage flux gj(x, t) per unit length, running perpendicularly 
between core and yoke. Hence, the magnitude of the main flux ®(x, rf) passing through core 
and yoke varies with the distance x from the heel piece end. 

If the current i, in the wire winding undergoes variations, we obtain through electro- 
magnetic induction secondary or eddy currents in the electrically conducting parts of the 
magnetic circuit. According to Lenz’ law the eddy currents give rise to opposing magnetic 
fluxes, which tend to slow down the rate of change of the total flux through a core section and 
to speed up the rate of change of the primary current. 

We shall show the following. If the geometric pattern of the eddy current distribution 
remains constant in time, it is possible to replace the eddy currents as to their influence on the 
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Fig. 1. Idealized relay with magnetic leakage. 
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primary current and on the magnetic fluxes by an equivalent winding with constant parameter 
values (resistance R,, inductance L,, mutual inductances M and M,,), carrying a current 
i, = the entire circulating secondary current. General expressions for the parameter values 
will be given. 


Hence, under the condition mentioned, we can for computational purposes replace 
a relay with eddy currents by a relay without eddy currents but with a delaying wire winding. 
Thus we can employ the theory of such a relay, given in Part 2. It should be stressed that the 
equivalent winding is a fictitious winding, which we can use formally in our calculations 
even if we may not be able to realize it physically. 


Units and symbols 
As in the earlier parts of this study the rationalized MKS-system of units is used through- 
out. 


In this part we denote 


quantities pertaining to the core by the superscript ’. 

54 "primary (magnetizing) winding by the subscript 1. 

eddy current system and the secondary (equivalent) winding by 
the subscript 2. 


As far as possible the symbols of Parts 1 to 3 are employed. See Fig. 1 and the list of 
symbols at the end of the paper. 
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Definitions of the secondary current and of the parameters of the equivalent winding 


We assume that the eddy currents in core and yoke run in planes perpendicular to the 
direction of x, and that their paths are symmetrical with respect to the core axis and to the 
symmetry planes of the yoke. This means (see Fig. 1): 

the eddy current paths in the core are concentric circles; the current density J; is a function 
of the radius @ of the path and of its distance x from the heel piece end, 


I, = 1,(@, x) (1 a) 
the eddy current paths in the yoke are parallel straight lines, with a direction parallel to 


the yoke surfaces; the current density J’, is a function of the distance y from the centre plane 
and of the distance x, 


Iz = 13 (y, x) (1b) 


Of course, J and J are in general also functions of the time ¢. 
Hence, we have a circulating current in the core, 
f doas (2a) 
x=0 


and in the yoke 


x=0 y=0 
We now define the current in the equivalent winding, the secondary current iz, by the con- 
dition that it should be equal to the entire circulating current, 
ip = ig + ig (3) 
We define the resistance of the equivalent winding, the secondary resistance Ry, by the 
condition that the heat developed when the current i, passes through R, be equal to the 
entire eddy current heat, 
112 
R,i2= + | 4, dt (4) 


core yoke 


We define the inductance of the equivalent winding, the secondary inductance L,, by the 
condition that the magnetic energy arising when the current i, passes through L, be equal to 
the entire magnetic energy of the eddy currents. Hence, introducing the magnetic vector 
potentials A, and Aj of the eddy current system, we shall have" 


2 [5 I,- A, dt +f; I, -Ag dt (5) 


core yoke 


i 
q 
q 


We finally define the mutual inductance M between the primary wire winding and the 
equivalent winding by the condition that the mutual magnetic energy of the primary current 
i, through the wire winding and the secondary current i, through the equivalent winding be 
equal to the mutual magnetic energy of the primary current and the eddy current system. 


Now, as is well known, the mutual energy can be expressed in two different ways, leading 
to two formally different expressions for M. 

If we use the vector potentials Aj, Aj of the primary current magnetic field and the eddy 
current vectors Ij, we have 


core yoke 
Here the integrals should be taken over the regions which carry eddy currents, i.e., over 
the volumes of the core and the yoke. 
If, instead, we use the vector potential A, of the eddy current magnetic field and the vector 
I, of the primary current, we have 


M ii. = (6 b) 
primary 
winding 
Here the integral should be taken over the region which carries the primary current, i.e., 
over the copper volume of the primary wire winding. 


Our use of the vector potentials in expressing the magnetic energies brings great advantages 
to our treatment inasmuch as we then only have to integrate over the current-carrying 
regions, while if we use expressions of the type { 1/2 B+ H dr, the integration must be extended 
to the entire volume of the magnetic field. 


Integration of (6b) 

It is easy to solve formally the integral in (6b). It yields the expression for M given in (9), 
which will prove to be of great importance to us. 

In carrying out the integration, we neglect the pitch of the winding, thus assuming that the 
current paths are circles, running perpendicularly to the core axis. We also assume that the 
primary winding has a very small thickness ¢ (see Fig. 2). This last assumption simplifies the 
calculations but is otherwise not essential. 
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Introducing in (6b) 


dt = 2mre dx 
el 
we obtain 
N, 
Mig= dx (7a) 
x=0 


where As, denotes the component of A, in the direction of I,, that is, in a direction tangential 
to the circumference of the core section. 

Now, an integral § A,ds being equal to the magnetic flux linked with the closed path of 
integration, we have 


(r, x) (8) 


where ¢’(r, x) = the total eddy current core flux at the distance x from the heel piece end. Thus, 
from (7a): 
1 
Mig= x) (7b) 


x=0 


Evidently, the last integral is equal to the total eddy current magnetic flux ®,, linked with 
the primary winding. Hence, 
M ig Dy (9) 
Thus, this formula, well-known from the theory of filiform circuits, is valid for our equi- — 


valent winding in the sense that the flux Mi, which the current in the equivalent winding sends 
through the primary winding, is equal to the eddy current flux ®,). 


Validity of the equivalent winding 

Evidently, the values of the quantities R,, L,, and M, which follow from (4), (5), and (6a) 
or (6b), are determined by the eddy current distribution at the moment in question. If our 
equivalent winding shall be practically useful, these values must be constant, independent 
of time. It is easy to see that, in accordance with our earlier prediction, R,, L,, and M are 
constant if the geometric distribution of the eddy currents is independent of time, i.e., if we can 
write 

x) =f(o, x) 12% 


with the same function y(t) in both expressions . 


If R,, L,, and M are constant, we can prove by energy considerations in conjunction with 
the relation (9) that the equivalent winding gives us the correct values of the currents i, and ig. 
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We start with the power relation 


dw 
Ey i, = + Reig) + (10) 
expressing that the battery power E£, i, is used up partly as Joulean heat R, i?+ R, i2, partly 


for increasing the magnetic energy 


1 1 


Employing (11), we can write (10) as 


di di di di 
Now, for the primary circuit we have the relation 
di d® 
+ ( - (13) 
or, making use of (9): 
di di 
+L,— +M—-E,=0 (14) 


dt dt 


The expression within the first pair of brackets in (12) being thus = 0, it follows that the 
expression within the second pair of brackets must also be = 0, i.e., 


Evidently, (14) and (15) are the equations governing the system formed by the primary 
circuit and the short-circuited equivalent winding. Hence, this coupled system yields correct 
values of the primary current i, and the total circulating eddy current i. 


The armature end flux 


We are above all interested in the total-armature end flux ®;. Following the procedure 
introduced in Part 2, Chapter 3, we express ®, by the aid of two mutual inductances, viz., 


M,, =the mutual inductance between the primary winding and a fictitious winding of 
one turn at the armature end, 


M,, = the mutual inductance between the equivalent winding and the fictitious armature 
end winding. 


dt dt 


®; is produced partly by the primary current i,, giving the contribution M,,i,, and partly 
by the eddy currents. Making use of the notation y’(r, x) of eq. (8), we can express the 
contribution from the eddy currents as 9’(r, 1). We write this 


1) (16) 


and take (16) as the definition of Mj. It is evident that M,, also is constant if the geometric 
distribution of the eddy currents is independent of time. 


We thus obtain for the armature end flux 
= + (17) 


i.e., (cf. Part 2, eq. (34)) the same formula as for a relay with a delaying wire winding. 


A physical picture of the equivalent winding 
Although it is hardly possible to materialize the equivalent winding, we can without 
difficulty form a physical picture of it. For simplicity we assume eddy currents in the core 
only. The general eddy current path then consists of a circle with the radius @ around the 
core axis. 
Fig. 3 shows two perpendicular sections of the core, one through the axis and one at right 
angles to it. 
We make use of the following symbols, dropping for the moment the superscript ’ and the 
index 2. 
V = core volume (volume element dr). 
S = surface, on one side of core axis, of a longitudinal core section (hatched in Fig. 3; sur- 
face element da). 
I =ccurrent density. 
U = 2xoE = voltage around an eddy current path (E = J/y = electric field strength). 
gy = 220A = total eddy current flux, linked with an eddy current path (A = vector potential 
of eddy current magnetic field). 
® = 2n0A, = primary current flux, linked with an eddy current path (A; = vector potential 
of primary current magnetic field). 
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According to the definition of i, we have 


[ (18) 


KY 


Further, eq. (4), defining R, gives us 


2 


or 
Rein = (19) 
Ig 
In the same manner we obtain from (5) and (6a): 
Ig 
Ss 
Mi, = f @ (21) 
2 


Now let us assume (see Fig. 3) that we divide the surface S into a large but finite number n 
of elements da, chosen so that the same current J da passes through all of them. We then obtain 
n circular turns with radii @ between 0 and r, all carrying the same current Ida = i,/n. 


We number the turns from 1 to n and denote by U;, 9,, and ®, the values of U, y, and 
®, pertaining to turn No. k. Then 


dp, d®,, 
and 
n d n d n 
2 
Egs. (19), (20), and (21) give us the approximate relations 
Rgig=- Leig=- (24a, b, c) 
ny ny 
Hence, from (23): 
diy di, 
Ryig +M (25a) 


i.e., the relation (15) already deduced. 
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Since the 7 turns all carry the same current i,/n, they can be connected in series. We then 
obtain a coil with a certain resistance R3, a certain inductance L3, and a certain mutual 


n 
inductance M*. This coil will have a voltage between its ends of > U; and it will link the 
1 


n n 
magnetic fluxes and > %,. Hence 
1 1 
Ra = M*i, (26a, b, c) 
1 1 


Comparing (24 a, b, c) and (26 a, b, c) we find 


R,=nR,; Lj=nL,; M*=nM (27a, b, c) 
From (23) follows 
tM (25 b) 


We can now summarize our investigation as follows: 


Pursuing the eddy current paths, we divide the entire core into a great number of turns 7, 
each turn carrying the same current i,/n. We insulate the turns from each other by very thin 
insulating layers, cut them open, connect them in series, and short-circuit the coil thus obtained. 
This coil will have a resistance R3, an inductance fod and a mutual inductance M*. For the 
short-circuited coil the relation (25 b) is valid. 


We could use this coil as our equivalent winding but for one serious disadvantage: in the 
limit n - co we would obtain a coil with infinitely large parameters, carrying an infinitely 


small current. 


Therefore, it is more convenient to employ an entirely formal equivalent winding with the 
finite parameters 


R* * M* 
R,= lim =; L,=lim =F M=lim — (28a, b, c) 


and carrying the finéte current i,. If we short-circuit this winding, eq. (25 a) = (15) holds. 
This equation differs from (25 b) only by the factor 1/n. 


The necessity of the condition that the geometric distribution of the eddy currents should 
be independent of time now becomes clear. This condition ensures that the m turns of our 
coil always carry equal currents, so that they can be connected in series, forming a fixed coil 
with constant parameters. Unfortunately this case is an exception; normally the eddy current 
pattern varies from one moment to the next, so that the configuration of the coil would have 
to change continuously. In other words, in most cases the equivalent winding cannot be 
used, because its parameters would vary with time. This is the reason why the equivalent 
winding is of no value in a study of the relay breaking problem. 
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On the other hand, the equivalent winding yields useful approximate results in such cases 
where the magnetic field produced by the eddy currents is small as compared with the main 
field produced by the current in a magnetizing wire winding. This latter field can then be 
considered as the only cause of the eddy currents, which can therefore be calculated from a 
homogeneously distributed magnetic flux varying with time. The parameters R,, Lo, and M 
will then be constant in time and easy to compute. 


If the eddy current magnetic field, although not small, has but little influence on the primary 
current and on the total magnetic flux, the same method can be reckoned upon to yield 
acceptable values of these two quantities. We shall make use of this when treating in the 
next chapter the relay closing problem. In this case we cannot expect very good results in the 
first moments where the eddy currents cause a quite large rate of increase of the primary 
current. 


CHAPTER 3 


Parameters of the eddy current equivalent winding 
of a relay on closing 


Distribution along the core of the primary current magnetic flux 


A constant primary current i, gives rise to a magnetic flux ®, in core and yoke, the magni- 
tude of which, on account of the magnetic leakage flux vy, varies with the distance x from the 
heel piece end. As shown in Part 1, the function ®, = ®,(x) can with sufficient accuracy 
be approximated with a parabola (Part 1, eqs. (28) and (20 a)). 


We introduce the quantities 


p 21 
=—: 
b (29a, b) 
+ a cosh 
(31a) 
+b(cosh  - 1) 


Putting 


12 


| 
i 
x 
=u 
= 


we can write the results of Part 1: 


D,(x) = B,(0) A(u) (32a) 
with 
A(u) =1 5 gtut (32b) 
and 
N, iy 
P,(0) = (32c) 


Generally 7?<1. Then 


(31b) 


1 
at+b+ab+- (1+ ) 
7? 

+56 67 (1 +55) 


grr 


In our subsequent deductions the parabolic function A(u) gives rise to the following easily 
calculated average values: 


m(u)= fA(u)du=u(1 + au/2 g?u?/6) (33a) 
n(u) =2 fm(u) da=w(1 + au/3 — g?u?/12) (33 b) 
m, = m(1) =1+4a/2-g?/6 (34a) 
n,=n(1) = 1 +a/3— 92/12 (34b) 
a= du =1+a-— (g?- a®)/3 ag?/4 +g¢4/20 (34d) 
n= fa +au)A(u)du=1 2a*)/6 ag?/8 (34e) 
0 
t= 3 = 1 + a— (7g?/4 a®)/5 — ag?/8 + 24/56 (34f) 


The eddy currents 


In calculating the eddy currents, we proceed as stated in Chapter 2, i.e., neglecting their 
own magnetic flux, we assume that they wholly result from the primary flux caused by the 
current i, through the magnetizing winding. This flux is uniformly distributed across core and 
yoke and varies parabolically along them; it changes with time in the same way as j,. 
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Hence, it can be written 


P, (x, 1) =A(u) 


where the function f(t) need not be specified. 


1) The eddy currents in the core. 
Applying Faraday’s law to a circular path of radius 9, we obtain 


2 


or, introducing ®,(x, ft) from (35) and putting — df(t)/dt = g(t): 


A(u) g(t) 


2ar 


For the entire circulating current in the core we obtain 


1 r 


x=0 e=0 


2) The eddy currents in the yoke. 


Here we apply Faraday’s law to a rectangular path of length a and width 2y: 


p(y, x)_ 
2a (x, t) 


Thus, 


BY. 


and the entire circulating current 


ef 
mic/2 , 
ig = dy y dx = g(t) 


a 
x=0 y=0 


3) The eddy current densities, expressed in the entire circulating current. 


(35) 


(36a) | 


(37a) 


(38a) 


(36b) 


(37b) 


(38b) 


For the entire circulating current—by definition also the current in the equivalent winding 


—we obtain 


ml c/2 
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We introduce the two quantities 


6’ and 6” show how the current i, is divided between core and yoke: 
6'+6"=1 (41a, b, c) 


Employing now (37 a, b), (39), and (40 a, b), we obtain the following final expressions for 
the eddy current densities: 


26’ 


The eddy current magnetic field 


The calculation of L,, M, and M,, requires that we know the magnetic field of the eddy 
currents. We need (see Fig. 4) the following two quantities: 


y’(o, x) = f vu’ H'(@, x) 220 do = the eddy current magnetic flux in the core, linked with a 
0 
circular eddy current path with the radius 9, situated at the distance x from the heel piece 
end. 
y 
y’(y, x) = {uv H"(y, x) 2a dy = the eddy current magnetic flux in the yoke, linked with a 
0 
rectangular eddy current path with the length a and the width 2y, situated at the distance x" 
from the heel piece end. 


heel piece air gap armature ' 
/ lpr) 
"| core | 
| 
YES — — Hy y|c 
| dx. Le -| 


Fig. 4. To the calculation of the eddy current magnetic fluxes. 
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We also introduce 


®,(x) = the total eddy current flux in core or yoke at the distance x from the heel piece end. 


Evidently 
PD, (x) =p (r, x) = @"(c/2, x) (43) 


In Part 1, Chapter 2, and again in Part 2, Chapter 3, are given the relations which govern 
the magnetic fluxes ® and ¢, produced by a constant current i, through a magnetizing wire 
winding covering the whole length of the core. These relations consist of two differential 
equations (Part 2, eqs. (62) and (63)) and two end conditions (Part 2, eqs. (64) and (65)). 


In our present case, where the fluxes are produced by the eddy currents, the first differential 
equation, expressing the flux continuity, can be taken over unchanged: 


(44) 


The same holds good for the end conditions: 
2d, (0) + s,pa(0) =0; 2, ®, (1) saga(l) =0 (45 a, b) 


In Part 2 we used the symbols ® and ¢;. In (44) and (45 a, b) we have replaced ® by ®, 
but for simplicity we have retained ¢y, although in the present case this flux also derives from 
the eddy currents. 


There remains the second differential equation, 


dp, Ny 
sP+ = (46) 


obtained as an application of Ampére’s law. This equation implies a uniformly distributed 
flux , caused by a current i, circulating on the core surface. In our present case this condition 
is not fulfilled, whence the equation cannot be taken over unchanged. 


In fact, the flux ®, is caused by eddy currents circulating within core and yoke, and the 
flux distribution across a section of core or yoke is strongly non-uniform. We shall show that 
(46) must be replaced by (51). 


We apply Ampére’s law to the rectangular path 1234 in Fig. 4, situated in the interspace . 
hetween core and yoke. No current being linked with the path, we obtain 
saga(x + dx) — saga(x) + {H(r, x) + H"(c/2, x)}dx=0 
or 


H'(r, x) + H"(c/2, x)= - (47) 


where H’(r, x) and H’(c/2, x) are the magnetic intensities at the surfaces of core and yoke. 
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We next apply Ampére’s law to the two rectangular paths enclosing the hatched areas 
within the core and the yoke, and obtain 


H(o, x)- H'(r, x)= i( > 0; x) (48 a) 

H"(y, x) - H"(e/2, x) =i"(>y, x) (48b) 

where i’( > 0, x) = the total circulating current per unit core length outside a cylinder with 
radius 

i’( > y, x) = the total circulating current per unit yoke length outside a rectangular 


strip with the transverse dimensions a and 2 y. 
Using (42 a, b), we immediately calculate 


We now turn to the two expressions for ®, given in (43): 
r c/2 
= x) 229 do= 3) 20 dy 


0 0 


Introducing here H’(o, x) and H’(y, x) from (48 a, b), and making use of (49 a, b), we find 


H'(r, x) =s' ®,(x) - A(u) ig (50a). 
H" (c/2, x) =s"®,(x) - A(u) - is (50b) 


Introducing these expressions in (47) and observing the relations s’+ s” =s, 6’+ 6”=1, 
we obtain the following equation, replacing (46): 


1 ” 
+0 


(51) 


Our next task is to find ®, = ®,(x) from (44), (45 a, b), and (51). Comparing (51) and 
(46), we see that for the purpose of computing ®, we can still assume a current on the core 
surface. This current, however, is not independent of x as in Parts 1 and 2 but varies para- 
bolically along the core. 

In principle, an exact integration of the system (44, 51) offers no difficulties. Still, we prefer 
an approximate method, which, while being sufficiently accurate, involves much simpler 
calculations. 
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The term s@®, in (51) represents the magnetic voltage drop per unit length of core + yoke, 
a quantity which is generally quite small. We can therefore, without substantial loss of 
accuracy, replace it by an approximate expression, which is known to us. With s®, known, 
an integration immediately gives us y, from (51), and by a subsequent integration we com- 
pute ®, (a better value than the one used in s®,!) from (44). Finally, the end conditions 
(45 a, b) give us the two integration constants. 


As our approximate value of s®, we choose the value, obtained by replacing the paraboli- 
cally distributed current in (51) with a constant current, giving the same total mmf. We are 
thus brought back to the relation (46) so that ®, follows at once from the results of Parts 1 


and 2. 
We denote the constant current per unit core length by ij. Hence, equalizing the mmf’s: 


2m 1 
0 
i.e. 
. 
lo 21 Ig (5 2) 
Substituting now in (32 c) ®,(0) for ®,(0) and i3/ for N,i,, we find 
S®, (0) = igl 
and hence, observing that S = /s: 
1+6°/3 
s®,(x) = s®,(0) A(u) = A(u) is (53) 
Using this expression for s®,, we obtain from (51): 
de 
=A(u)-i3 (54) 
with 
ae mn? \ 
= = 5 
‘3 2m! (: g* +7? ‘2 
Integrating (54), while observing (33 a), we find 
| 
ga(x) = pa(0) m(u) (56) 


We now integrate (44), making use of (56), (33 b), and the relation s; = /S,. It follows that 


= x93 (0)- m(u) (57) 
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The unknown constants ®,(0) and ¢,(0) are determined by (45 a, b). One finds that (57) 
can be written 


| 


®,(x) = ®,(0) (1 + au) (58) 
with 
®,(0) = (59a, b) 


Knowing ®,, we can now determine the magnetic field intensities H’(o, x) and H’(y, x). 
The eqs. (48 a, b) together with (49 a, b) and (50 a, b) give us 


H(o, x) =s' B, (x) -2 (2) (60a) 
H"(y, 3) + - ‘baw: (60 b) 


Introducing ®,(x) from (58), and paying attention to (55) and (59 a, b), we finally obtain 
(1 +5) (1 [ea n(a)| + 
2m 1 3 +n? 2 
+8 (2) (61a) 
H"(y, x)= ( +5) ( wt) + 


From the expressions (61 a, b) we compute without difficulty the fluxes y’(o, x) and 
y’(y, x). The result is 


| 


1 +6"/3 mn? y 


B(i +au)-= | + 


(=) | (62b) 
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As a check we compute ®,(x) from both expressions: 


®,(x)=9'(r, x)= 9" (c/2, x) = 


and observe that the same value of ®,(x) follows from (58), combined with (55) and (59 a, b). 


The resistance of the equivalent winding 
From eq. (4) follows 


i 
x=0 e=0 4 x=0 y=-c/2 ? 
Taking J,, 1 from (42 a, b) and performing the integrations, we easily find 
Q, [2% 16a/3c 
R,=—-(— 6°+—— 
( 8 (65) 


The inductance of the equivalent winding 


According to eq. (5) we have 
Lj (66) 


, 1 , , , ” 1 ” ” ” 
2 


2 
yoke 


with 


core 
Reasoning as when deducing (20), we obtain 


ig ig i, is 


x=0 


With the aid of (42 a, b) and (62 a, b) this gives us 


{(1+ 5) (1 (361-31) +05 (69a) 
6” 4 6” 1 


and, if we make use of 6’ + 6” = 1: 


4 6”\2 m,n? 1 62 166” 
{(1 +3) (1- Br, sh) +4 (= (70) 
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The mutual inductance between the equivalent winding and the primary winding 
For the mutual inductance M we derive two different expressions, one from (6 a) and the 


other from (6 b). 
Taking first (6 a), we proceed as when calculating L, and find 


M = M’ + M’ (71) 
with 

ig iy ip iy 

x=0 e=0 x=0 y=0 
Here 


®, = ®;(e, x) =the primary current magnetic flux in the core, linked with a circular eddy 
current path of radius @ situated at the distance x from the heel piece end. 


®{ = P{(y, x) =the primary current magnetic flux in the yoke, linked with a rectangular 
eddy current path of length a and width 2y situated at the distance x from 


the heel piece end. 
Introducing the total primary current flux 


(x) = (r, x) = Pi (C/2, x) 


we have 
x)=(—) Pi); Pi Y= P, (x) (73a, b) 
with ®,(x) according to (32 a, b, c). 
We now integrate (72 a, b). The result is 
= 2m, Si + 2m, Si (74 a, b) 
and 
6” 1 
75 


We pass on to (6 b). This relation has already been transformed to (7 b), which can be 


written 
1 


ig 
x=0 


Introducing y’(r, x) from (63), we find 


m= (145) (1 1+; (77) 
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Generally, (75) and (77) lead to somewhat different values of M. This is due to the fact that 
the expression (63) for the eddy current flux y’(r, x) depends on the approximate value of 
s®, which we used in (51). For normal values of S = sl, however, (75) and (77) yield very 
nearly the same numerical result. This confirms that our approximation of s®, was permis- 
sible. 


It is not difficult to see from (51) that in two cases, for S = 0 and for S, = oo, our method 
of solution is exact. By comparing (75) and (77) one also finds that in both of these cases the 
two expressions for M coincide. : 


The mutual inductance between the equivalent winding and a fictitious winding of one 
turn at the armature end 


Eq. (16) immediately gives us 


My, =" (78) 
2 
where ¢’(r, 1) follows from (63). It ensues 
6” mn? 1 
(1 +5) (1-25) (1 (79) 


CHAPTER 4 


The closing of relays without magnetic leakage and with eddy 
currents in the core only 


In this chapter we study a simplified case, accessible to an exact treatment. It allows us to 
judge the usefulness of the equivalent winding solution by comparing it with the exact 
solution. 


The exact solution 


a) Expansion theorem solution for large t. 


As mentioned in Part 3, a paper by I. Heruirz®, published in 1920, implicitly contains the 
solution of the eddy current problem for a relay with eddy currents in the core only, and 
without magnetic leakage. Making use of Heaviside’s expansion theorem, Herlitz could 
express his quantities by means of exponential function series valid for not too small values of 
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time ¢. His expression for the magnetic flux, written with the symbols of our telephone relay 
problem, is presented in Part 3, eq. (29). 
We give below the expansion theorem solutions for the following four quantities: 
i, = i,(f) = the primary current. 
® = P(t) = the magnetic flux. 
H = H(t, e) = the magnetic intensity in the core. 
I = I(t, e) = the eddy current density in the core. 


(@ = distance from the core axis.) 


The subsequent notations are used: 


and 
A;= A(x) == 
2 J, (xi) 


where J, and J, are Bessel functions and x; (i = 1, 2, 3, ...) are the (real and positive) 
roots of 


A(x) - 2 4.K=0 
R,S't 
The solutions are: 
() 
i,(t E, 
= = oj} > — 80 
i=] 
80b 
i=1 
H(t, 0) (1+K)A x) Nii 
»0 “a i H. =H 100 80c 
i=l 
I(t, @) A;+K ) y 
= IL=I 0, 80d 
h > ©. on 


i=1 
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b) Series solutions for small t. 

For small values of t we employ series solutions in powers of 1/2, obtained by developing 
the operators in powers of p-/*, The calculations are similar to those given in Part 3 for a 
relay with magnetic leakage (cf. Part 3, eq. (124), for the armature end flux; with leakage the 
series is in powers of f1/4), 


The following expressions for the primary current i, and the magnetic flux ® contain terms 
up to and including ¢* in i, and ®/t: 
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va\r 2 4 
1/2 1/2 3/2 
D, 3Va\r 8 t 8 


In these formule 


and 
A,= (k-4) 
A,= 2gk +3/8 
A, = — (k® — 3gk? + +3/4)k— 3/8) 
Ay=  (k*- + 3/8)k?- (¢ +1)3k/4 + 63/128 
A,= — (k® — Sqk4 +(q? + 1/4) — + 9g/4 + 9/8)k? + (q + 3/2) 3k/4 - 27/32 


with 


The equivalent winding parameters 


Our present, simplified problem is characterized by the conditions 


Si = 0; y”=0 
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They give rise to the following numerical values: 


=m =h=1 


Introducing the total magnetic reluctance 


4+S8°=2+2,4+5 


and observing the relations 
Ss. B 1 


+7? Sj, 00 Sy aot 2; 


we find from (65), (70), (75) or (77), and (79) the following values of the equivalent winding 


parameters: 
2x 1 3S” 
82a, 
N,. 1 
(82c, d) 


We add the parameters of the primary winding: 


2 
(82e, f) 


The formule for primary current and armature end flux (Part 2, eqs. (44) and (46); see 
also the list of formule at the end of this paper) make use of the following quantities: 


primary time constant 


L, N2/R 
(83a) 
| secondary time constant 
| L | 38’ 


time constant in the formula for the armature end flux 


Mo, M yl 
leakage factor 
2 
(83d) 
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Basic data of the quantitative study 


The quantitative study of this and the next chapter refers to the relay already investigated 
in the preceding Parts 1, 2, and 3. The data of this relay were given in Part 3, Chapter 4 (for 
the geometric data see the second paragraph following the chapter heading; for the electric 
and magnetic data see Tables 2 and 6). 


In short, the data refer to an Ericsson RAB relay with a Lancashire iron yoke and a Lanca- 
shire or silicon iron core, an air gap of 0.10 or 0.87 mm, and a winding of 10000 turns with 
the resistance 560 Q. In the present chapter we put the leakage reluctances S$, = oo. 


Two series of permeability values are employed, viz., the series of high values, 
corresponding to the mmf 300 At, already used in Part 3, and a series of low values, equal to 
1/10 of the values of the first series. 


In Table 1 the different cases treated are specified and suitably designated. 


Table 1. The numerical cases treated in Chapters 4 and 5. 


Air gap Relative 
Case Core mm permeability 
L87(u) Lancashire 0.87 5 400 
L10(u) 0.10 3 400 
$87(u) Silicon 0.87 4 500 
$10(u) 0.10 3 400 
L87(u/10) Lancashire 0.87 540 
L10(u/10) 0.10 340 
$87(u/10) Silicon 0.87 450 
$10(u/10) 0.10 340 


Results 


Fig. 5 gives primary current and magnetic flux vs. time in a case where the eddy current 
influence stands out clearly. The curves for actual relays with magnetic leakage have the same 
general appearance. 


Without eddy currents both curves are exponential with the time constant 13.7 ms. With 
eddy currents the flux curve starts with an unchanged rate of increase. As compared to the 
exponential curve, it then slows down, remaining all the time below this curve. In other 
words: the eddy currents always cause a decrease of the flux value. 


The current curve starts with an infinite rate of increase; we find from the power series 
development that di,/dt ~ const. x t-"/? for very small values of t. The curve at first rises 
well above the exponential curve but after a while its rate of increase becomes so small that 
it cuts across the exponential curve, remaining subsequently below it. 


Thus, the primary current curve displays a pronounced knee. In Part 2 we found the same 
phenomenon in relays with delaying slugs and, in its extreme, in the sleeved relays, where 
the current rises to its final value almost at once. 


26 


t 
| 
| 
| 
| 
\ 
| | 


no eddy 


oat 


fe) 10 20 30 40 50 60 70 80 ms 
t 


X 11199 


Fig. 5. Primary current and magnetic flux on closing a relay without 

magnetic leakage and with eddy currents in the core only. Solid curves: 

case L87(u/10), Lancashire core, relative permeability of core material 

= 540, air gap 0.87 mm. Dotted curve: exponential curve, time constant 
13.7 ms, valid for relay without eddy currents. 


The effects of a change in the permeability value are of interest. An increased permeability 
of course leads to a larger time constant and hence to an exponential curve with a slower 
increase. Closer study (cf. Figs. 10 a to 13 a) shows that an increase of the permeability 
entails a decrease of the eddy current influence and a shift toward higher time values of the 
intersection between the current curve and the exponential curve. 


The four pairs of figures, Figs. 6 a, b to 9 a, b, give an idea of the accuracy of the equivalent 
winding solutions. They all refer to the large air gap (0.87 mm), the one of primary interest 
at closing. The results for the small air gap (0.10 mm) are similar but somewhat more 
favorable as to the applicability of the equivalent winding. 


Figs. 6 a, b refer to case L87(), the case with the most marked magnetic skin effect and 
the largest error in the equivalent winding solution. For four time values, 2.5, 10, 40, and 
160 ms, Fig. 6 a gives the distribution across the core section of the magnetic intensity H (to 
the left), and the eddy current density / (to the right). 


One sees that still at 160 ms the distribution of the magnetic flux across the core is 
extremely inhomogeneous and the eddy current density quite far from the linear increase 
from core axis to core surface required by our theory (cf. eq. (42 a)). Yet, as shown by 
Fig. 6 b, the equivalent winding solution gives reasonably accurate results except for the 
primary current at small times. 

At the other end we have case $87(u/10), Figs. 9 a, b, with a fairly small magnetic skin 


effect. Here the eddy current density follows the linear law much better, and the equivalent 
winding gives acceptable results down to quite small time values. 
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Basic data of the quantitative study 


The quantitative study of this and the next chapter refers to the relay already investigated 
in the preceding Parts 1, 2, and 3. The data of this relay were given in Part 3, Chapter 4 (for 
the geometric data see the second paragraph following the chapter heading; for the electric 
and magnetic data see Tables 2 and 6). 

In short, the data refer to an Ericsson RAB relay with a Lancashire iron yoke and a Lanca- 
shire or silicon iron core, an air gap of 0.10 or 0.87 mm, and a winding of 10000 turns with 
the resistance 560 Q. In the present chapter we put the leakage reluctances S,; = 0. 


Two series of permeability values are employed, viz., the series of high values, 
corresponding to the mmf 300 At, already used in Part 3, and a series of low values, equal to 
1/10 of the values of the first series. 


In Table I the different cases treated are specified and suitably designated. 


Table 1. The numerical cases treated in Chapters 4 and 5. 


Air gap Relative 
Case Core mm permeability 
L87(u) Lancashire 0.87 5 400 
L10(u) 0.10 3 400 
S87(u) Silicon 0.87 4 500 
$10(u) 0.10 3 400 
L87(u/10) Lancashire 0.87 540 
L10(u/10) 0.10 340 
$87(u/10) Silicon 0.87 450 
$10(u/10) 0.10 340 


Results 

Fig. 5 gives primary current and magnetic flux vs. time in a case where the eddy current 
influence stands out clearly. The curves for actual relays with magnetic leakage have the same 
general appearance. 

Without eddy currents both curves are exponential with the time constant 13.7 ms. With 
eddy currents the flux curve starts with an unchanged rate of increase. As compared to the 
exponential curve, it then slows down, remaining all the time below this curve. In other 
words: the eddy currents always cause a decrease of the flux value. 


The current curve starts with an infinite rate of increase; we find from the power series 
development that di,/dt ~ const. x t-"? for very small values of t. The curve at first rises 
well above the exponential curve but after a while its rate of increase becomes so small that 
it cuts across the exponential curve, remaining subsequently below it. 


Thus, the primary current curve displays a pronounced knee. In Part 2 we found the same 
phenomenon in relays with delaying slugs and, in its extreme, in the sleeved relays, where 
the current rises to its final value almost at once. 
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Fig. 5. Primary current and magnetic flux on closing a relay without 

magnetic leakage and with eddy currents in the core only. Solid curves: 

case L87(u/10), Lancashire core, relative permeability of core material 

= 540, air gap 0.87 mm. Dotted curve: exponential curve, time constant 
13.7 ms, valid for relay without eddy currents. 


The effects of a change in the permeability value are of interest. An increased permeability 
of course leads to a larger time constant and hence to an exponential curve with a slower 
increase. Closer study (cf. Figs. 10 a to 13 a) shows that an increase of the permeability 
entails a decrease of the eddy current influence and a shift toward higher time values of the 
intersection between the current curve and the exponential curve. 


The four pairs of figures, Figs. 6 a, b to 9 a, b, give an idea of the accuracy of the equivalent 
winding solutions. They all refer to the large air gap (0.87 mm), the one of primary interest 
at closing. The results for the small air gap (0.10 mm) are similar but somewhat more 
favorable as to the applicability of the equivalent winding. 

Figs. 6 a, b refer to case L87(), the case with the most marked magnetic skin effect and 
the largest error in the equivalent winding solution. For four time values, 2.5, 10, 40, and 
160 ms, Fig. 6 a gives the distribution across the core section of the magnetic intensity H (to 
the left), and the eddy current density J (to the right). 


One sees that still at 160 ms the distribution of the magnetic flux across the core is 
extremely inhomogeneous and the eddy current density quite far from the linear increase 
from core axis to core surface required by our theory (cf. eq. (42 a)). Yet, as shown by 
Fig. 6 b, the equivalent winding solution gives reasonably accurate results except for the 
primary current at small times. 

At the other end we have case $87(/10), Figs. 9 a, b, with a fairly small magnetic skin 
effect. Here the eddy current density follows the linear law much better, and the equivalent 
winding gives acceptable results down to quite small time values. 
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Fig. 6a. Distribution across the core section 

of magnetic intensity and eddy current density 

on closing. No magnetic leakage. Eddy currents 
in core only. Case L87(u). 
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Fig. 6 b. Comparison between the exact 
and the equivalent winding solutions for 
primary current and magnetic flux on 
closing. Deviations from values of relay 
without eddy currents, in percent of 
final steady state value. No magnetic 
leakage. Eddy currents in core only. 
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Fig. 7 a. Distribution across the core section 

of magnetic intensity and eddy current density 

on closing. No magnetic leakage. Eddy currents 
in core only. Case L87(u/10). 
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Fig. 7 b. Comparison between the exact 
and the equivalent winding solutions for 
primary current and magnetic flux on 
closing. Deviations from values of relay 
without eddy currents, in percent of 
final steady state value. No magnetic 
leakage. Eddy currents in core only. 
Case L87(u/10). 
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Fig. 8 a. Distribution across the core section 

of magnetic intensity and eddy current density 

on closing. No magnetic leakage. Eddy currents 
in core only. Case S87(). 
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Fig. 9 a. Distribution across the core section 

of magnetic intensity and eddy current density 

on closing. No magnetic leakage. Eddy currents 
in core only. Case $87 (u/10). 
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Fig. 8 b. Comparison between the exact 
and the equivalent winding solutions for 
primary current and magnetic flux on 
closing. Deviations from values of relay 
without eddy currents, in percent of 
final steady state value. No magnetic 
leakage. Eddy currents in core only. 
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Fig. 9 b. Comparison between the exact 
and the equivalent winding solutions for 
primary current and magnetic flux on 
closing. Deviations from values of relay 
without eddy currents, in percent of 
final steady state value. No magnetic 
leakage. Eddy currents in core only 
Case $87(u/10). 
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CHAPTER 5 


The closing of relays with magnetic leakage and with eddy 
currents in core and yoke 


Numerical results 


We now turn to the study of actual relays with magnetic leakage and with eddy currents 
in both core and yoke. We employ the values of the leakage reluctance S,, given i Part 3, 
Chapter 4, Table 2. As in Part 3 (see Part 3, Chapter 4, Table 6 and the two paragraphs 
above the table), we take account of the eddy currents in the armature and the heel piece, 
as well as in those parts of the core which are not covered by the winding, by properly in- 


creasing the values of the core and yoke conductivities. 
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Fig. 10. Influence of eddy currents 
on primary current and armature 
end magnetic flux on closing. 
Deviations from values in relay 
without eddy currents, in percent 
of final steady state value. Relay 
with magnetic leakage. 
Cases L87(u) and S87(w). 


Ai; % = CaseL 87(ul0) 
Thin lines:no eddy 


5 current 
in yoke 


X 11205 


Fig. 11. Influence of eddy currents 
on primary current and armature 
end magnetic flux on closing. 
Deviations from values in relay 
without eddy currents, in percent 
of final steady state value. Relay 
with magnetic leakage. 
Cases L87(u/10) and S87(u/10). 
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Still, the problem remains so complicated that it is hardly possible to deduce a usefu\ 
exact solution, that is, to find from the operational expressions given in Part 3 numerically 
useful time functions. We therefore make use of the equivalent winding. We can do this 
with great confidence. In fact, for the simplified problem of Chapter 4, the equivalent winding 
gave us satisfactory results down to at least 5 ms, and neither the magnetic leakage nor the 
eddy currents in the yoke can be expected to bring about any basic change of our theory. 


We have calculated, for the eight cases given in Table 1 of the preceding chapter, the 
equivalent winding parameters, and from these the armature end flux and the primary 
current on closing. For comparison the cases of eddy currents in the core only and of no 
eddy currents have also been worked out. 

The results are collected in the four Figs. 10 to 13. The curves are drawn on a logarithmic 
time scale. In each case they present the deviation of the actual current or flux value from 
the corresponding value in a relay without eddy currents. The deviations are given in percent 
of the final steady state value. 
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Fig. 12. Influence of eddy currents 
on primary current and armature 
end magnetic flux on _ closing. 
Deviations from values in relay 
without eddy currents, in percent 
of final steady state value. Relay 
with magnetic leakage. 
Cases L10(u) and S10(u). 


Fig. 13. Influence of eddy currents 
on primary current and armature 
end magnetic flux on closing. 
Deviations from values in relay 
without eddy currents, in percent 
of final steady state value. Relay 
with magnetic leakage. 
Cases L10(u/10) and $10(u/10). 
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Some general conclusions 


The eddy current influence on current and magnetic flux in actual relays on closing is quite 
small. Therefore, it is not easy to compare our calculations with experimentally deter- 
mined oscillograph records. From such a comparison it appears, however, that of the two 
series of permeabilities in Table 1, only the one with high values leads to values of current 
and flux which are close to the actual ones. By a suitable adjustment of these high values, 
one could no doubt obtain a very satisfactory agreement with the oscillograms. 

In other words, during the closing operation the low initial permeability is of little im- 
portance. One should not find this astonishing. To be sure, for small time values the magnetic 
flux is small but its distribution across the core and yoke sections is strongly inhomogeneous. 
The flux density and consequently the permeability will therefore reach large values in those 
parts of the iron which effectively transmit the flux and where the permeability value is 
hence of importance. One can presume that the permeability value corresponding to the final 
steady state should yield good results down to fairly small time values. 

In our study of the breaking operation in Part 3 the same argument was put forth, the 
flux on breaking being small but highly inhomogeneous for high time values. 

We can therefore state that both on closing and on breaking, a constant permeability value, 
in the neighbourhood of the (generally high) value corresponding to the D.C. steady state of the 
relay, will allow a good representation of the entire course of events. In both cases the perme- 
ability values for small magnetic intensities are of little interest! 

It appears that the facts discussed here have not always been duly understood. In the 
literature on telephone relays a homogeneous flux distribution is often assumed, whereas— 
on closing as well as on breaking—the magnetic flux in core and yoke is in general strongly 
inhomogeneous during the major part of the transient period. It should be observed, 
however, that when the flux passes over to the air gap, it almost immediately assumes 
a distribution which is very close to the one characterizing the static case. Hence, the 
inhomogeneity in the iron can be assumed to have only little influence on the pulling force 
on the armature. Actual force measurements support this conclusion. Its elucidation by 
further theoretical study would be of value. : 
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Principal symbols and formule 


t = time. 

1 = length of core and yoke (length of primary winding). 

x = distance along core and yoke, counted from the heel piece end. 
r = radius of core. 

a = width of yoke. 

c = thickness of yoke. 


’ = electric conductivity of core material. 


y” ata ” ” ” yoke ” 

yp’ == magnetic permeability of core material. 

” ” ” yoke ” (In Part Y25 My, Hy.) 


N, = number of turns of primary winding. 


S’ =s'l = reluctance of core. 

=] = », yoke. (In part 3: S,, Sy.) 
S =S’ + S’ = sl = reluctance of core + yoke. 
2) = reluctance of heel piece end. 

4, = Sy 5, armature + air gap. 

S, = s,/l = total leakage reluctance. 

R, = resistance of magnetizing (primary) circuit. 
L, = inductance ,, 9 9 ” 

R, = resistance of equivalent (secondary) winding. 
L, = inductance ,, 


M = mutual inductance between primary and equivalent winding. 
M,, = mutual inductance between primary winding and a fictitious winding of one turn at 


the armature end. 
M,2, = mutual inductance between equivalent winding and a fictitious winding of one turn 


at the armature end. 
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E, =constant emf in magnetizing circuit. 

i, = = current in primary winding. 

ijoo = £,/R, = steady state value of primary current. 

ig = i,(t) = current in equivalent winding (total circulating eddy current). 
@ = (x, t) = total magnetic flux in core and yoke. 

®,(t) = (1, t) = magnetic flux at the armature end (in Parts 2 and 3: ®,). 
®,, = M,, E,/R, = steady state value of armature end flux. 

Pi = 9i(x, t) = magnetic leakage flux per unit core length. 


Formule 


Auxiliary quantities 


2o, 2, 

S;’ S;’ = Si 


ab +acoshyn +b 


+b (cosh 7 1) 


1 1 
a+b+ab+ +56) 
—b+—n? 


m,=1+a/2—- g?/6 

ny, = 1 +a/3 - g?/12 

Py =1+a/4- g?/20 

1 a2)/3 ag*/4 + 24/20 
r,=1+a-(g?—- 2a?)/6 ag?/8 


t, =1+a- (7g?/4 - a®)/5 — ag?/8 + 24/56 


1 
m, +5 bn, 


a+b+ab 
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(29a, b) 


(30a, b, c) 


(31 a) 


(31b) 


(34a) 
(34b) 
(34c) 
(34d) 
(34e) 
(34f) 


(40a, b) 


(59b) 


7) 


Equivalent winding parameters 


(i (75), cf. eq. (77) 
My, = (1 +5) (1 [ ec +a)- (79) 


For S;= 00, y” =0 (no magnetic leakage, eddy currents in core only), see eqs. (82 a—f), 
(83 a—d). 


Time constants (cf. Part 2) 


L L M,, M 
T,=—; Ty=T-— 
V(2 +T,\?_ of, T, 
where 
Primary current and armature end flux (cf. Part 2) 
i t,- "-(t,.- ™ E, 
100 T— Te 
t t 
(%,- Tye Tne E, 
—=1 ®,, = D;(~) = M,, 
oo ( ) 11 R, 
Gothenburg, February 1961 
Institute for Theoretical Electricity and Electrical Measurements. 
Chalmers University of Technology. 
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